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We discuss some properties of zero and near-zero modes of the Dirac operator, as observed in a 
recent simulation of 2-flavor QCD. The quarks have been implemented with the so-called Chirally 
Improved Dirac operator, which obeys the Ginsparg-Wilson relation to a good approximation. We 
present geometrical and statistical properties of these eigenmodes and eigenvalues. 


1 Motivation and introduction 

Throughout the last decades of lattice QCD topological excitations have played a crucial 
and sometimes controversial role. It is natural, in particular in view of the purpose of this 

volume in honor of Adriano Di Giacomo, that they are a topic also hereP 

The topological charge of a gauge configuration is a well-defined concept only in the 
continuum and for differentiable fields. In quanfum field fheory fhe gauge configuralions 
are in general non-differenliable. Laffice discrefizafion furns ouf fo be an ideal fool for 
doing a non-perfurbafive gauge invarianf regularizafion of QCD; however infroducing a 
smallesf disfance fhus we lose fhe continuum in such sfudies. In facl, a consfrucfive def¬ 
inition of fhe fOTological charge is only possible when cerfain smoofhness requiremenfs 

are introduced^ If is also quite demanding to implement in actual calculations. 

For that reason various other methods have been addressed, among them ultra-local 

definitions of the topological charge density together with various modifications® 
When considering the total topological charge of gauge field configurations cooling tech¬ 
niques have been employed in the attempt to smoothen the gauge configuration and get 
rid of discretization artefacts like dislocations for example. 

One of the cleanest tools for identifying the topological charge is through the measure¬ 
ment of a costly fermionic observable: eigenmodes of the Dirac operator. The Atiyah- 

Singer index theorem® relates the topological charge of a gauge configuration with the 
difference between the number of left-handed and right-handed zero modes. Unfortu¬ 
nately the discretization in itself complicates this approach. In the continuum zero modes 
also have definite chirality. Dirac operators on the lattice violate explicitly chiral sym¬ 
metry, as it is defined in fhe continuum. Only recenfly one has rediscovered a relation® 
generalizing fhe continuum chiral symmefry fo fhe laffice by adding a local piece, fhaf 
vanishes in fhe continuum. In its simplest form this so-called Ginsparg-Wilson (GW) 
relation reads 

{D,D^ = DD^, ( 1 ) 

and it may be associated with a symmetry transformatiorH, that becomes the chiral trans¬ 
formation in the continuum limit. This symmetry protects the zero modes, which now 
have definite chirality even on the lattice. 

Dirac operators obeying the GW relation may thus be utilized to identify the topo¬ 
logical charge of a gauge configuration by identifying the number and chirality of the 
zero modes. Various GW type Dirac operators have been introduced: Domain Wall 


FermionJ^I^ISin] the Perfect FermionJ^l^or the Chirally Improved (Cl) FermionsO^I^ 
They all fulfill the GW condition not exactly but in some approximation. They are also 
more complicated to deal with than the simple Wilson operator, having many more cou¬ 
pling terms or introducing an extra di mension . The only known explicit and exact re¬ 
alization, the so-called overlap operator^I^, is even more expensive to implement in 
computer simulations. 


With better chiral behavior there comes an associated extra advantage. The simple 
Wilson Dirac operator has scattered real modes near the origin in the complex plane. 
These become true zero modes in the continuum limit. At present day lattice spacing 
values they lead to spurious zeros even at non-zero mass of the quark fields (see, however. 

Ref. ^1^. Improving the chirality also improves this situation. GW type fermions have 
the chance to allow for a better approach to the chiral limit, i.e., to come closer to the 
physical value of the pion mass without having to go to very large and fine lattices. Exact 
GW fermions have no spurious modes, i.e., no eigenvalues below the values of the mass 
parameter. 


Real zero modes may thus be used to obtain the total topological change; its fluctua¬ 
tion gives the topo logical susceptibility which, in the quenched case is related to the mass 

of the rj' (cf. Refs .^^^^^^. For the dynamical case its dependence on and m has been 
discussed recentl}^222E2l2l] 


So-called near zero modes are associated with complex conjugate pairs of small eigen¬ 
values. Spontaneous chiral symmetry breaking relates their density to the value of the 

chiral condensate via the Banks-Casher relationP^l Random matrix theory (RMT) de¬ 
scribes the probability density of such eigenvalues for operators belonging to various 
universality classes including that of the QCD Dirac operator. Su ch a pro bability density 

depends on the topological sector, the quark mass and the volumeP^J^ 


For all these reasons zero modes are of high interest: both, leading to important phys¬ 
ical observables, and as a technical tool to identify the behavior and quality of the simu¬ 
lation. 


Table 1: Parameters for the simulations; the first column denotes the run, for later reference. The gauge 
coupling is /3i. the bare quark mass parameter am, HMC-time denotes the length of the run (number of 
trajectories), and the lattice spacing has been determined via the Sommer parameter. 


# 

X T 

Pi 

am 

HMC-time 

05 [fm] 

a 

12^ X 24 

5.2 

0.02 

463 

0.115(6) 

b 

12^ X 24 

5.2 

0.03 

363 

0.125(6) 

c 

12^ X 24 

5.3 

0.04 

438 

0.120(4) 

d 

12^ X 24 

5.3 

0.05 

302 

0.129(1) 

e 

8^ X 16 

5.3 

0.05 

1245 

0.135(3) 























2 Dynamical fermion simulation 


The Cl Dirac operator Dqi was constructed by writing a general ansatz for the Dirac 
operator 

16 

Dij = Y.a\^{U)Tk, ( 2 ) 

k=l 


where (A: = 1... 16) are the 16 elements of the Clifford algebra and a^j{U) are 
linear combinations of path ordered products of links U connecting lattice site i with site 
j. Inserting D into the GW relation and solving the resulting algebraic equations for 
the coefficients of the linear combinations yields Dqi- In principle this can be an exact 
solution, but that would require an infinite number of terms. In practice the number of 

terms is finite and the operator is a truncated series solution to the GW relationl^^J 
Cl fermions have been already extensively tested in quenched calculations (see, e.g.. 


Ref In these tests it was found that one can go to quark masses below 300 MeV 
without running into the problem of exceptional configurations (spurious zero modes). 
On quenched configurations pion masses down to 280 MeV could be reached on lattices 
of size 16^ X 32 ( lattice spacing 0.148 fm) and about 340 MeV on 12^ x 24 lattices. 

In recent workl^^l^ we have studied the Cl fermions in a dynamical simulation of 

Eol 


QCD with two light flavors. All technicalities are discussed in Ref.l^. For definiteness 
we just mention that our gauge action is the Liischer-Weisz actionl^, that we used stout 


smearing^of the gauge fields as part of the Dirac operator definition, and that the Hybrid 
Monte Carlo method was implemented to deal with the dynamics of the fermions. The 
lattices were (up to now) of moderate size: 8^ X 16 and 12^ x 24, for lattice spacings 
between 0.11 and 0.14 fm. Table [D summarizes the simulation parameters of the runs 

nol 

discussed here (see, however, ReTi^for a more complete list). 


3 Zero modes and near zero modes 

3.1 Isosurface plots 

When dealing with zero modes some visualization of the geometric structure may be de¬ 
sirable. The instantons, which provide finite energy analytic solutions to the Yang-Mills 
field equations, are prime candidates for the topological excitations. In some scenarios 

a gas of instantonsl^ supposedly leads to chiral symmetry breaking. In this picture the 
low-lying near-zero modes come from overlapping instanton - anti-instanton pairs. 

A possible analysis tool is to introduce the gauge invariant density of the correspond¬ 
ing eigenvectors ipcaix), in particular 

Po{x) = ^f>ca{x)fca{x) , P5{x) = 'fca{x){'y5)apfcp{x) , (3) 

c,« c,q,/3 

with the normalization J2x Po{x) = 1 (c denotes the color index and a,/3 the Dirac 
indices). The sum J2x P^i^) gives the chirality (V'ItsIV') of tho mode. For exactly zero 
modes of an exact GW operator one has po{x) = 3zp^{x) locally. 
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Figure 1: Iso-surfaces at time slices close to the maximum of the scalar density. Rows 1 and 2: typical 
case of a single lump; row 1: po (I = 15.6), row 2: ps (7 = 12.1). Rows 3 and 4: several lumps in one 
eigenvector; row 3: po (I = 3.2), row 4: ps (7 = 2.4). 

For approximate GW operators like ours, the real eigenvalues of D are not exactly 
zero, but the corresponding eigenmodes are the only ones with non-vanishing chirality. 
As a typical example let us study two eigenvectors corresponding to slightly different 
but real eigenvalues of one configuration. In Fig. ^ we compare iso-surfaces at po = 
0.25 of the maxima values, shown as 3-cuts at time slices close to the density maxima. 
Comparing the density values for pQ with ps we find that the chirality of the mode is 
indeed located mainly where the scalar density is concentrated. The two sets of plots 
also show typical distributions: the zero modes are as often located in individual blobs 
as they are distributed on several such lumps. Note, that we discuss here eigenvectors for 
individual real (“zero ”) modes and not sums of such. 

The inverse participation ratio 

h = VY.Pk{xf (4) 

X 

(V is the lattice volume in units of the lattice spacing) gives a global information on the 
localization of the eigenmode. Its value ranges between V denoting point localization 
and 1 for an uniformly spread mode. If a mode is uniformly distributed on N sites it 
has I = V/N, thus l/I indicates the fraction of the volume occupied by the mode. 
Fig.|2lshows the histogram for the inverse participation ratio measured on the zero mode 
eigenvectors for the run (a). The values of Iq for the zero modes are quite similar to the 

results for quenched simulations for comparable lattice size and spacingP^The average 
values for run (a) are Iq = 14.0(3.4) and Is = 8.7(2.0). 

In Ref.it was observed that the lumps may change their positions when the bound¬ 
ary conditions change. A possible interpretation is that the lumps carry fractional topo¬ 
logical charges. We normally use anti-periodic boundary conditions in time direction. 























Figure 2: Histogram for the inverse participation ratios Iq and Jsfor the zero modes in run sequence (a). 


However for a subset of 10 of our configurations we also determined the eigensystem for 
periodic b.c.; in one of those we observed a change in the number of zero modes. For the 
other configurations the number stayed the same and the density showed only very little 
change. 

Tunneling between different topological sectors appears to be a problem for HMC 
implementations of the overlap action; various intricate methods have been suggested to 
deal with it. We do not seem to have such a problem and observe frequent tunneling. 
We determined the eigenvalues and thus the topological charge only for every 5th con¬ 
figuration. Comparing the number of changes of the topological sector we get (for the 
runs (a,b,d) and (e)) between 12 and 19 such changes along a HMC-time distance of 100, 
without obvious correlation to the run parameters. These values are a lower bound for the 
actual number of tunneling events. As is well-known from experience with other Dirac 
operators, when evaluated over longer HMC-periods than t he one s available to us, the 

tunneling frequency may show much longer correlation timeP^^^ One ma y ex pect then 
longer wavelength fluctuations in the topological charge than observed here!“ 


3.2 Density of near zero modes 

Random matrix theory gives the density distributions of the low-lying non-zero eigen¬ 
values in universality classes depending on the general symmetry properties of the Dirac 
operator. The densities should scale with a scaling variable proportional to V Im(A). The 

proportionality constant may be related to the fermion mass and the chiral condensateP^^^ 
Although we have limited statistics, in particular for the larger lattice, we still identify 
such volume scaling as can be seen in Fig. |5| 

In order to exhibit the scaling behavior we plot the abscissa scaled according to the 
respective volume. We also indicate the positions of the mean values of the histograms 
and find good agreemenf wifh fhe expecfed volume scaling wifhin fhe limifed sfafisfics. 


3.3 More on the smallest eigenmodes 

Recenfly fhe disfribufion of fhe smallesf eigenvalue of fhe hermitian Wilson Dirac op- 

erafor has been sfudied in a dynamical fermion simulafionESl for large laffices and fhe 
sfandard Wilson Dirac operafor. There if was found fhaf fhe median of fhaf disfribufion 
is near fhe value of fhe so-called AWI-mass, fhe physical quark mass defermined from 
fhe PC AC-relation. The widfh of fhe disfribufion appears fo scale proportional fo ajsfV 
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Figure 3: Density of the Im(A) values of the 1st, 2nd and 3rd lowest non-zero modes for the runs (d) and 
(e) at /3i = 5.3 and am = 0.05 for two volumes (top: V 12 = 12® x 24, bottom: Vs = 8® x 16) in the 
sectors of topological charge i/ = 0 and 1 / = 1. The abscissa of histograms are scaled with the volume. 


(in physical units). This width gives information on the range of mass values accessible 
in such simulations. The eigenvalues occurring to the left of the AWI-mass lead to the 
so-called spurious zero modes, i.e., a zero of the Dirac operator at positive quark masses 
and thus a spurious singularity in the quark propagator. 

Let us briefly consider the situation for a perfect GW operator, in particular one fol¬ 
lowing the simple form Eq. ^ for zero mass. The eigenvalue in that case lie on a circle 
with radius 1 and center at 1 in the complex plane. One may add a mass through 

/ jYl \ 

Dcw{m) = ( 1 “ y j Dgw + ml . (5) 

The smallest (in terms of distance from the origin) eigenvalue of Dcwim) is either real 
Amin = m or on the circle at, say Amin = m ± ie. “ The probability distribution of the 
density of |Amin| at finite volume consists of a delta-function at |Amin| = m (from the 
exact zero m odes) and the folded distribution of the lowest non-zero mode, known from 

rMtP^ 


“The corresponding smallest eigenvalue of the hermitian GW operator 'ysDawi'm) is 














































Figure 4: Sketch of the distribution of the smallest 
eigenvalue |A| for a GW Dirac operator D(yif^{rn) 
on a finite lattice. The vertical thick line at A = m 
symbolizes the delta function due to the real (“zero”) 
modes. The shaded area represents the modes on the 
circle, which dominate in the infinite volume limit; 
in this limit, the width of this part of the distribution 
shrinks to zero oc 1/V^. 
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Figure 5: Distribution of the absolute value of the smallest eigenmodes for run (a) (l.h.s.) and (b) (r.h.s). 
The vertical, dotted lines indicate the values of the AWI-mass as determined in Ref. The quark mass has 
been converted to physical units using the lattice spacing determined from the Sommer parameter. 


Fig. m gives a sketch of the situation for an exact GW operator. The weight of the 
delta-function is proportional to the number of configurations with at least one exact zero 
mode and thus related to the finite volume topological susceptibility. Assuming the topo¬ 
logical charge |i/| has Gaussian distribution with (zz^) = xV the probability of the zz = 0 
sector is P = 1/\/27r(zz^) cx; 1/ \/V ; this defines fhe normalizafion of fhe nearesf non¬ 
zero modes disfribufion. (Nof e, however, fhaf fhe suscepfibilify for dynamical fermions 

vanishes wifh fhe fermion The weighl of fhe della-funclion is 1—P. 

In Fig. |5l we show fhe disfribufion obfained for fwo of our simulafion runs. Only 
fhe eigenvalues below fhe value of fhe AWI-mass are spurious modes. The sifuafion is 
similar fo fhaf of an exacf GW operator. In fhe sfudy for dynamical Wilson fermions the 
median of the corresponding distribution was close to the value of the AWI mass. In our 
case, the AWI-mass is much below that median. This motivates the use of GW type Dirac 
operators. 


4 Summary 

We have implemented the Chirally Improved Dirac operator, obeying the GW condi¬ 
tion to a good approximation, in a simulation with two species of mass-degenerate light 
quarks. Studying the low-lying eigenmodes allows us to check the quality of the Dirac 
operator and the simulation. Here we surveyed some details concerning the geometrical 
properties of these eigenmodes and the probability distribution. 
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